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ABSTRACT 

Standard magnetohydro dynamic theories, such as the mean field dynamo theory, have 
been criticized when the back reaction of the magnetic field on turbulent motions is ne- 
glected. For the dynamo, this back reaction has been argued to suppress the turbulent 
motions required for optimal mean field production. Here it is suggested that if the mag- 
netic field is spatially intermittent, for example residing in flux tubes, the back reaction 
on turbulent flows may be significantly reduced. 
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1. Introduction- The Mean Field Dynamo (MFD), an elegant theoretical mechanism 
that allows a large scale magnetic field to grow exponentially at the expense of shear and 
small scale turbulent energy, has been studied to explain the origin of magnetic fields in 
astrophysical bodies such as planets, stars, and galaxies [1,2,3,4]. The difficulties with 
dynamo physics highlight fundamental issues of magneto-hydrodynamics that are still not 
well understood. In particular, the effect of magnetic fields on turbulent flows. 

In principle, the MFD growth of a large scale magnetic field in a differentially rotating 
system occurs as turbulent motions induce formation of magnetic loops from a seed field 
[1,2,3]. As a result of the coriolis and buoyancy forces, the turbulence is cyclonic; the loops 
in both hemispheres twist in the same direction, creating a large scale loop of poloidal 
field. The turbulent diffusion of the outer portions of these mean field loops ensures that 
the net flux of mean field lying in the region of interest is non- vanishing [1]. Differential 
rotation shears the large scale poloidal field, generating a large scale toroidal field. The 
new field then incurs the same loop forming process, providing the feedback for exponential 
growth. The strength of the field is limited by the available turbulent energy. Numerical 
simulations of "kinematic" dynamos in which the back reaction of the field growth on the 
turbulent eddies is neglected, can produce magnetic topologies consistent with observations 
of stellar and galactic magnetic fields [5,6]. 

In reality, the small scale, root mean squared (rms) field energy density grows much 
more rapidly than the mean field [1], violating the kinematic approximation. Though, for 
example, observations of the solar photosphere and the dispersed heavy element distri- 
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bution in the Galaxy indicate the presence of reasonably uninhibited turbulent motions 
in amidst equipartition magnetic fields [1], theoretical studies [7,8] suggest that the back 
reaction should suppress these motions and thereby inhibit the dynamo. As first discussed 
by Piddington [9], the argument proceeds as follows: As an eddy tries to displace a parcel 
of fluid threaded by a field loop in equipartition with the turbulent energy density, the 
loop acts as a rubber band, restricting further motion, and recoiling the parcel to its ori- 
gin. Statistically, transport of material is significantly reduced. Simulations support this 
intuition by finding that an increasingly large fraction of plasma motions are locked into 
oscillating modes rather than zero frequency (diffusive) modes as the magnitude of the 
initial ordered field is increased [8]. 

Here it is suggested that turbulent motion may survive the back reaction if the mag- 
netic field is concentrated in flux tubes. Because the tubes' Alfven speeds can be larger 
than the eddy velocities, reconnection between tubes can be rapid. This would reduce the 
back reaction force on the turbulent velocities. The diffusion and helicity coefficients of 
the MFD equation would be reduced from the kinematic theory only by the fraction of a 
typical tube thickness which does not reconnect in an eddy turnover time. (In the special 
case for which the magnetic field is totally composed of topologically unlinked loops or 
cells, interchange motions between the cells would not be restricted by magnetic forces 
and reconnection would not be required for diffusive motions.) 

First, the derivation of the dynamo helicity and diffusion coefficients is outlined and 
the specific role of the kinematic approximation is highlighted. The stages of dynamo 
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growth are discussed and the formation and role of flux tubes is then addressed. An 
estimate of the dynamo coefficients is given based on the physical ideas presented. Finally, 
a similar role of intermittency for star forming regions is mentioned. 

2. Dynamo Coefficients and Interpretation of Approximations- The magnetic field 
and velocity are taken to be B = b + B and V = v + V respectively, where b and v are 
fluctuating quantities with zero mean and B and V are mean quantities. The Reynolds 
relations [11] are also required. These are dt^(KiHj) = (dt^(KiHj)), and (Kihj) = 0, 
where Ki = Ki + ki and Hj = Hj = hj are components of vector functions of position x 
and time t, the brackets indicate the mean value, and <9t iX is the derivative with respect 
to x or t. For ensemble averages, these relations hold when the turbulence is correlated 
on time scales short relative to the variation time scales of the mean quantities. For the 
spatial average, defined by (K(x, t)) = \(\~ s f**£ K(s, t)d 3 s, these relations hold when the 
average is taken over a large enough scale, that is when / << \(\ « L, where L ~ B/VB 
is the scale of the mean field variation, and / ~ b/Vb ~ v/Vv. 

The induction equation derived from the non-relativistic Maxwell equations is [1] 

8(B + h)/dt = Vx [(V + v)x (B + b)] + z/ M V 2 (B + b), (1) 

where vm is the magnetic viscosity. Astrophysical magnetic Reynolds numbers are large 
and the last term in (1) is unimportant on the energy containing eddy scales for destroying 
magnetic field energy (but it does provide dissipation on the smaller scales allowing a 
turbulent cascade, and is important locally, at the intersections between thin flux tubes.) 
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Taking the average of (1), ignoring the last term, gives the MFD equation 

dB/dt = Vx[Vx B + (vxb)]. (2) 

The turbulent electromotive force (TEMF) given by (v x b), is written [2] 

(v x b) = aijBj + p ijk VjB k , (3) 

where is the helicity dynamo coefficient, and flijk is the turbulent diffusion dynamo 
coefficient. Under the assumption of isotropic turbulence oiy = dija and flijk = eijh$- 

Working in a frame for which V = 0, subtracting (2) from (1) gives 

8b/dt = Vx(vxB) + Vx[vxb- (vxb)] + b- VV. (4) 

Usually, for simplicity, the 3rd and 4th terms of (4) are neglected straight away, which 
is a procedure called the "first order smoothing (FOS) [1,2]" approximation. Instead, I 
will keep these terms and then show what this approximation means. Plugging (4) into 
(3), the dynamo coefficients become expansions of time-ordered exponential series [13] in 
powers of r c /r e d, where r c is the turbulence correlation time and r e d ~ l/v is the eddy 
turnover time. Under the assumption of isotropy, taking the first terms in the series give 
the standard forms [1,2] 

a= (-l/3)(v(t)- / dt'Vxv(t')), (5) 
Jo 

P=(l/3)(v(t). /'dfv(f)>. (6) 

Non- vanishing a means non- vanishing helicity, essential for mean field growth as described 
above. That (3 acts as a diffusion coefficient for the mean field is evident from plugging (6) 
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into the TEMF in (2). For spatially homogeneous turbulence, (3 becomes the coefficient of 
the diffusion term on the right side. 

Note that dropping the higher order terms that led to (5) and (6) is essentially the 
equivalent of the FOSA. The inclusion of higher order terms is an additional complication 
which is extraneous to, and independent of the focus on the back reaction in this paper: 
The back reaction reflects the specific effect of the magnetic field on the velocity, not the 
combination of velocities that appear in the dynamo coefficients. 

In the usual kinematic approximation, v is prescribed independently of b and B. This 
motivates the use of (4) to eliminate b in (3). However, the functional forms of the dynamo 
coefficients are exactly valid even when v = v(b, B). The specific v(b, B) will depend on 
the application, but the field would always be inhibitive to turbulent motions. Note that 
because the time scale for growth of the small scale field is much shorter than that of the 
mean field [1] , the most important back reaction comes from the small scale field. This is 
the natural interpretation of studies [7] which show effects of the back reaction for values 
of the mean field much less than equipartition with the turbulence. 

3. Phases of Dynamo Growth and Flux Tube Formation- Observations of the sun 
[14,15] and simulations of MHD turbulence in low f3 ave (= P pa rt/Pmag where P par t and 
Pmag are the average particle and magnetic pressures) plasmas [16] indicate that the mag- 
netic field tends to concentrate in flux tubes. Intermittency in magnetic field strength in 
the Galactic interstellar medium is seen as well [17]. Determining the size of flux tubes 
and the role of f3 ave will be addressed later but how such intermittent structures form is 
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addressed first. 



A working mean field dynamo would incur phases given an initial seed field [13]. In 
the 1st phase, turbulent energy stirs up the rms random magnetic field to equipartition. 
In the 2nd phase, in principle, the mean field also nears equipartition with the turbulence 
and/or shear while the small scale field remains at or near equipartition. In the 3rd 
phase, the dynamo works to sustain the mean field. The back reaction is straightforwardly 
unimportant only during the 1st phase which lasts for a time ~ r e d- 

Flux tubes can form in phase 1 and Ref. [18] is relevant. There, the evolution of a 
seed magnetic field in homogeneous, isotropic turbulence is studied. An important result 
of [18] is that the field tends to concentrate locally into flux tubes or ropes. The field 
only grows in a local region when turbulence conspires to produce to proper stretching, 
twisting and folding [19]. Assuming that the mean field is constant over an eddy turnover 
time scale, d'B/dt ~ db/dt so that the induction equation for the total field can be used 
to explore the growth of the random field for phase 1. In Lagrangian form 



As in [1], consider the initial Bi(0) to be aligned with a line segment Sxi(0). Then at all 
later times Sxi(t) is aligned with Bi(t). Eq. (8) then implies that the length of a line 
segment parallel to the field satisfies 



where f(x,t) is a random function of the turbulent velocities. Although (l(x,t)) can be 
shown to increase exponentially [1], l{x,t) is equally likely to decrease or increase at a 



dBi/dt = BiViVj. 



(7) 



dl(x,t)/dt = f(x,t)l(x,t), 



(8) 
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position x. Similarly, the field would be equally likely to increase or decrease there, and 
thus a natural spatial intermittency would result. The rms field in this picture, grows 
nearly to equipartition with the volume averaged turbulent flow energy density, which 
likely happens in only of order a time ~ r e d [1] for a modest seed field. The back reaction 
is most inhibitive after equipartition ensues, and this is therefore the case of interest for 
this paper. 

The growth of field in a particular localized region results from stretching and dragging 
of seed flux by a turbulent flow. As the bundle or tube is dragged favorably for exponential 
growth, material will be inhibited from seeping into the tube from the ends since there the 
field is weak and the force on particles opposes the direction of field line convergence [20] . 
The amount of mass in a tube should thus remain the same, or decrease. Assuming that 
it remains the same, incompressibility implies a constant volume. Thus the cross sectional 
area of the tube -kt\ oc Z -1 , and from flux freezing B t oc /. The edge of the tube is a current 
sheet, with no flow normal to it, and thus pressure is balanced across it : 

Pt + l = (p ext + l)(B ext /B t ) 2 , (9) 

where (3 t = 87rP partjt /B^ and (3 ext = 87rP partjext /B^ xt with B ext is of order the initial seed 
field, P P art,t is the tube particle pressure, and P par t,ext is the external particle pressure. 
Eq. (9) shows that (3 t « (3 ext since B t » B ext in equilibrium. 

Since the tube pressure is balanced by the external pressure, the magnetic energy 
density can be higher than the turbulent energy density at the tube locations when (3 t is 
small, and then the volume filling fraction of the tubes would necessarily be small. To 
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see this note that the average magnetic and particle pressures satisfy P mag ~ fPmag,t and 
Ppart ~ (1 — f)P P art,ext where / is the fraction of volume occupied by magnetic flux tubes, 
Pmag,t is the magnetic pressure in the flux tubes and P pa rt,ext is the particle pressure 
external to the tubes. Thus (3 ave ~ [(1 - f)/f]P pa rt,ext/Pmag,t = (1 + 0t)(l - /)//. Thus 

//(i-/) = (i + &)//W 

Each energy containing (outer) scale eddy of wavelength / stretches a tube to length / 
and radius rt- The thickness of each tube, rt, can be estimated by balancing the magnetic 
and turbulent eddy drag forces [9,21]. This gives 

(B?/4nr c )(nr*) ~ C dPext vf2r t , (10) 

where B t is the magnitude of the field in the flux tube, p ex t is the density outside the flux 
tube, and Cd is the coefficient of turbulent drag. Since / is a wavelength, the radius of 
curvature r c can be estimated by 1/4 when the tube maximally responds to the turbulence. 
In equilibrium, B^/8n ~ P pa rt,ext(^ + Pt), so that (10) gives, for (1 - /) ~ 1, 

r t /l = CaYMf/A-K ~ 0.06M Z 2 (1 + t ), (11) 

where V is the adiabatic index and Mf = vf /(TP ext /p ext ) ~ (3~^ e ~ f(l + I3 t )~ l , when 
equipartition between turbulent and magnetic energy is assumed. For the last similarity in 
(11), Ca was estimated from the "drag" crisis [21] which reduces Ca < 1 at large turbulent 
Reynolds number Ri. Assuming Ri > 1000, Cd ~ 0.4. 

4- Role of Flux Tubes- Once dynamo growth enters the 2nd and 3rd phases, the back 
reaction of the field on the turbulent eddy becomes important with respect to the net 
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transport of magnetic field. In particular, the rapid growth in magnetic tension inhibits 
it from traveling much more than ~ /. Unless reconnection with another tube can happen 
before the tension response, the tube will react back. Both the diffusion and helicity 
coefficients require inhibited motions of the turbulent velocity. For example, inhibited 
turbulent diffusion would mean that an ink mark on some tube statistically incurs zero 
net displacement from its initial location (i.e. oscillating motions) instead of increasing its 
separation from the initial point with time. 

No inhibition of motions would mean that a tube could reconnect with a partner in 
the time it takes to move a distance rt, namely a time rtjv\. After a reconnection, the tube 
would change one of its ends, and would then move in a different random direction from 
whence it came, before impacting another tube. The process would continue, enabling 
for example, an effective diffusion. The reconnection time scale is r t F(R,M,t)/vA,t where 
F(RM,t) is the function of the magnetic Reynolds number associated with the length scale 
rt, and VA,t is the tube Alfven speed. There would be no inhibition to turbulent motions 
when 

v A , t / Vl ~ M-\l + A)" 172 ~ 0111/(1 + (3 t ) 1 / 2 > F(R M , t ), (12) 

where the second similarity follows from equipartition. For Sweet-Parker (SP) reconnec- 
1 li 

tion, F(R,M,t) ~ Rm t wnue for Petschek (PK) reconnection F(RM,t) ~ Lni?M,t [1]- Note 
that R M ,t ~ (r t /l)R M ,i = OmR M ,iM^(l + ft) ~ 0.06# M ,z(l + Pt)/Pave from (11), where 

Rm,i is the standard magnetic Reynolds number associated with /. For the more strin- 

l li 

gent SP case, (12) then becomes f3 ave /(l + (3 t ) > 0.25-R^ t which is likely satisfied in or 
above the solar convection zone of the sun [1]. Note that if the field were diffuse and not 
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concentrated in flux tubes, then then (3 ave ~ j3 t and the inequality in (12) could not be 
satisfied. 

Note that (11), and the line following it, imply that r t /L < f . But each tube fills a 
fraction ~ r 2 / L 2 < f 2 of an eddy volume. Thus there are > 1// tubes per eddy volume 
when (1 — /) ~ 1. 

5. Application to the Dynamo Coefficients- The simplest way to describe the effect 
of flux tubes, is to say that fast reconnection significantly reduces the back reaction terms 
of the Lorentz force on the velocity field in the equation of motion. Thus the velocities in 
the dynamo coefficients would be approach their kinematic values the more efficient the 
reconnection. In the presence of equipartition magnetic flux tubes, the reduction from their 
kinematic limit would be determined by the fraction of tube Lorentz force that contributes 
to the back reaction, namely the fraction that cannot reconnect during r e d- Equivalently, 
using (13), and noting that dynamo coefficients depend on two powers of the velocity, we 
have 

OLhr ~ Min \5tkini 

and 

for ~ Min kin , p ave (l + Pty^kin/FiRMj) 2 ], (14) 

where the subscripts kin and br refer to the kinematic values and the values including 
the back reaction, respectively. The right sides of (13) and (14) are the minima of the 
quantities in brackets. 

It can be useful to think of the diffusion coefficient for, as measuring the fraction of 
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eddy energy per mass, contained in motions which random walk rather than oscillate. For 
a given amount of total eddy energy, a stronger back reaction means a higher fraction 
of non-zero frequency modes [8]. To see this, note that when the velocity is given by 
a stationary random field, (v(t) ■ v(t + r)) = C(r). Then the Fourier transform gives 
C(ui) = f_ oo d,TExp[iwt](v(t) ■ v(t + r)). The zero frequency component, by comparison 
with (6) then satisfies f3 = (l/4)(7(0). where C(0) is the Fourier transform of the velocity 
correlation (i.e. the power spectrum of the velocity field) evaluated at zero frequency. The 
amount of eddy energy per mass in the zero frequency modes (i.e. the fraction contributing 
to the diffusion coefficient) is the non-zero contribution to the diffusion coefficient 

6. Discussion Rapid reconnection, resulting from a concetration of magnetic fields 
into low fit regions, may overcome the back reaction on turbulent motions. Diffusion of the 
mean field and helicity would be enabled not necessarily by removing a large amount of field 
energy on the outer scale, but by allowing individual flux tubes to avoid recoiling back to 
their points of origin. The most important motions would be enabled on the energetically 
dominant scales. However, a steady forcing of the turbulence on these outer scales would 
give a cascade to small scales as in (e.g.) [22], maintaining a constant magnetic + turbulent 
energy density on the outer scale. The irreversible dissipation would occur on the smallest 
scales. Eqns (13) and (14) apply most effectively when f3 ave » Max[l, (3 t \. The value of j3 t 
determines how effectively flux tubes evacuate and there may be a non-linear dependence 
on Pave- Actual values will have to await future simulations. 

If the shear energy were much larger than the turbulent energy and could be dumped 
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into the field fast enough, the magnetic energy may exceed the turbulent energy. Then the 
first similarity in in the line following (11) would become < and the inequality condition 
between the last terms in (12) would be stricter than required, since the third term would 
be < the first two. 

Finally, note that a similar rapid reconnection between evacuated tubes in low (3 ave 
star forming molecular cloud regions of the ISM may remove material from field lines and 
initiate collapse [23]. In an initially uniform f3 ave « 1 plasma, non-linear compressional 
Alfven waves clump material on the scale of the energy containing eddies to density en- 
hancements of order l/(3 ave . The Alfven speed associated with the sparse regions is large, 
allowing rapid reconnection, closed loops formation, and dissipation. Thus intermittency 
can also lead to fast reconnection in plasmas with (3 ave « 1. 

Thanks to G. Field for discussion. 
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